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ABSTRACT 

Approximately half of the matter in the Universe is “unbound” at z = 0, according to 
N-body simulations such as the Millennium Run. Here, we use the milli-Millennium 
simulation to examine the distribution of unbound matter in relation to the dark 
matter halos which host galaxies. We measure the unbound matter within two types 
of windows, using a halo dependent radius and a fixed radius at several different scales. 
We also consider the timescales over which a halo can accrete the local unbound matter 
at z = 2 and z = 0. Finally, we compare the unbound matter to observable properties 
of galaxies, such as local galaxy count environment and stellar mass. We find that 
halos at z = 2 can accrete far more of the nearby unbound matter over a Hubble time 
than halos at z = 0 and that 78% of particles within 5 Rvir of a halo at z = 2 will 
be accreted by z = 0, compared to 36% of particles within 5 h~^ Mpc of the halo. 
We also find that galaxy count environment is closely related to the amount of nearby 
unbound matter when measured on the same scale. 

Key words: galaxies: clusters, galaxies: evolution, galaxies: haloes, methods: statis¬ 
tical 


1 INTRODUCTION 

Galaxies are bright and easy to observe, but on their own 
they are not a complete picture of the Universe. They live 
in the densest peaks in the cosmic web and are surrounded 
and connected by diffuse matter. Estimates based on hydro- 
dynamic simulations predict that equal numbers of baryons 
at z = 0 are found in galaxies, warm-hot intergalactic mat¬ 
ter, and cold filamentary gas ( Dave et al.||1999 1. This diffuse 
gas has been observed at z ~ 2 — 3 in Ly-a emission ( |Can-| 
[talupo et al. 2014 Martin et al. 20141, at numerous red- 
shifts in Ly-Q absorption ( Lynds|1971 as well as in X-rays 
||Werner et al.|[2008| |Ma et al.||2009[ [Planck Collaboration| 
et al.||20T3l at z = 0 — 0.5. Radio observations will be pos¬ 


sible with future telescopes such as the Square Kilometre 
Array ( [Takeuchi et al.|2014[ ). 

The halos in which galaxies live are part of the cosmic 
web, a structure that has formed over the age of the universe. 
Structure formation in a ACDM universe is well defined both 
analytically ( Press &; Schechter||1974[) and through N-body 
simulations (e.g. Springel et al. 20051. Halo growth is as¬ 
sumed to happen through mergers and accretion events (e.g. 
|Lacey Cole|1993[|Wechsler et al.|2002p , with the division 
between the two being somewhat arbitrary. 

“Major mergers” are usually defined as having a mass 
ratio between the progenitor galaxies of Mi/M 2 < 4 (e.g. 
Lotz et al.pOllI although sometimes the mass ratio cutoff 


ratios greater than this cutoff are considered “minor merg¬ 
ers”, although sometimes these have an upper mass ratio 
limit of 10. For every decade increase in mass ratio, the 
number of mergers increases by a factor of 5, although this 
is limited to only a few decades by the mass resolution of 
the N-body simulation used ( Fakhouri fc Ma|2008 |. 

Very small mergers, or accretion evens, happen much 
more frequently than major mergers and contribute a sig- 
nihcant amount of mass to an evolving halo. Unresolved 
mergers can account for 40% of the total bound mass in 


a halo (Genel et al. 20101 although their relative impor¬ 


tance in comparison to major and minor mergers is some¬ 
what halo mass dependent ( [Fakhouri fc Ma|2010p . Likewise, 
Wang et al. ( 2011| find that the fraction of nearby substruc¬ 


ture that is bound to a halo depends mostly on halo mass, 
although the fraction of nearby substructure that is unbound 
to a halo is dependent on the local tidal field. Structure out¬ 
side the extent of the halo is also necessary to accurately 
measure the matter power spectrum ([van Daalen fc Schayej 


20151. 


Gas traces these additions, which can be very smooth 


(Angulo & White 20101, either undergoing a shock and heat¬ 


ing up to the virial temperature, or radiating most of its en¬ 
ergy away and accreting as a cold flow ( jBirnboim fc Dekelj 
2003 Keres et al.||2005 (. 


As Fakhouri & Ma (20081 discuss, N-body simulations 


can be as large as 10 (e.g. Mihos fc Hernquist||1994 1. Mass limit the mass ratios that can be explored by the resolu- 
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tion of the simulation. Halos still grow through accreting 
the unbound particles individually or as sub-resolution ha¬ 
los. Gas traces these additions, either undergoing a shock 
and heating up to the virial temperature, or radiating most 


of its energy away and accreting as a cold flow (Birnboim & 
Dekel|2003 Keres et al.||2005l. 


It is this pre-accretion unbound matter we are most 
interested in. A common question addressed by galaxy en¬ 
vironment studies is how does a galaxy’s environment affect 
its evolution. We turn this question around to study the 
opposite: how does a galaxy’s evolution affects its environ¬ 
ment. One method of measuring environment, which has 
been used extensively in the literature, is to count the num¬ 
ber of nearby galaxies within an aperture of fixed size. In this 
paper, we will continue to use this definition and explore its 
connection to the amount of unbound matter nearby. While 
both definitions can be classified as a galaxy’s environment, 
we will use the word “environment” to refer to the galaxy 
count within a fixed aperture. We will focus on three ques¬ 
tions regarding this unbound environment: 

(i) How is unbound matter distributed in the universe? 
And how does that change with time? 

(ii) How likely is the nearby unbound matter to accrete 
onto a galaxy? 

(hi) Are there correlations between the amount of un¬ 
bound matter around a galaxy and observable properties of 
that galaxy that we can use to make predictions? 

This paper is organised in the following way: in Sec¬ 
tion we discuss the simulation and model used. We break 
our analysis into three parts to address each of the ques¬ 
tions above. In Section we discuss the first question of 
the distribution of matter in the universe and its evolution 
with time. In Section]^ we look at the free-fall times of the 
unbound matter and determine how likely it is that matter 
will accrete onto a halo. In Sectionwe consider observable 
properties of galaxies and determine if there is a “smoking 
gun” of future potential accretion. In Section we extend 
our discussion of the three questions posed above. Finally, 
in Section[^ we summarise our findings. To avoid confusion 
between the gas in the intergalactic medium (IGM) and the 
dark matter it traces, we use the term diffuse to refer to 
the former and unbound to refer to the latter. We also use 
the term IGM to refer to all gas that is not bound to a 
halo, not differentiating between the phases as is common 
in hydrodynamic simulations (e.g. Dave et al.|2010 l. 


2.1 The Simulation 

For our N-body simulation, we use the milli-Millennium 
Run (mMR), a test ca se for the full Millennium Simulation 
(Springel et al. 20051. The mMR is a 62.5 h~^ Mpc co¬ 


moving periodic box with (270)® particles of mass 8.6 x 10® 
h~^ Mq. It was run with cosmological parameters from the 
WMAPl survey ( Spergel et al.|2003 |. The simulation tracks 
perturbations from 2 = 127 until 2 = 0 over 64 snapshots. 
We focus on the 32 snapshots from 2 = 2 until 2 = 0, look¬ 
ing mostly at these two redshifts specifically. The mMR is 
ideal for this work because, unlike in the full Millennium, 
the particle data is publicly available. 

We use the online database of halos and unbound par¬ 


ticles, as determined by the halo finder SUBFIND (Springel 
et al. [2001 1. SUBFIND is a friends-of-friends algorithm that 
requires a certain number of particles to be linked together 
for a halo to be considered bound. In this paper, we use the 
standard minimum number of 20 particles (giving a mini¬ 
mum halo mass of 1.7 x 10^® h~^ M©) and a linking length 
of b = 0.2. According to SUBFIND, 23% of the particles are 
considered bound at 2 = 2, and 49% at 2 = 0. Using AHF 
| |Knollmann fc Knebe|[2009| ), an adaptive mesh halo finder, 
38% of the particles are bound at 2 = 0. With Rockstar 


(Behroozi et al. 20131, a friends-of-friends algorithm with 


an adaptive linking length, 55% of the particles are bound 
by 2 = 0. When using 100 particles as the minimum halo 
mass at 2 = 0, the same three halo finders show 41%, 34%, 
and 48% for SUBFIND, AHF, and Rockstar, respectively. 

In other tests of N-body simulations with varying box 
and particle size, we find the fraction of bound particles 
to be far more dependent on the halo finder and chosen pa¬ 
rameters than on the simulation resolution, regardless of the 
resolution. This is due to the shape of the halo mass func¬ 
tion, which when integrated between lx Mparticie and 20x 
Mparticie (1-6. the unbouud portion of the distribution) is typ¬ 
ically 40-60% of the total mass in the box at 2 = 0 (at least 
for the resolution ranges cosmological simulations are cur¬ 
rent run at). It is a technical distinction and the reason why 
we separate the “unbound mass” (a simulation term) and 
“diffuse mass” (more physical). In the real Universe much 
of the diffuse mass is probably bound to their local struc¬ 
tures out to scales larger than Rvir. But in the simulations, 
as used by the semi-analytic models, this mass is thrown 
away. 

Halos are discovered as groupings of more than 20 parti¬ 
cles, but are re-defined in the following way. The virial mass 
(Mvir) and the virial radius (Rvir) are related by 


2 METHODS 

Halos, by construction, are definable objects with a finite 
extent in dark matter simulations. Unbound matter, how¬ 
ever, is a nebulous term that can refer to anything from sub¬ 
resolution halos to the dark matter traced by the intergalac¬ 
tic medium. For the purpose of this paper, sub-resolution 
halos and IGM-traced structure are considered to be the 
same. The unbound matter, therefore, consists of the dark 
matter particles that are not currently part of a halo. Like¬ 
wise, halos are easy to compare because they have specific 
properties, such as mass and extent. Unbound matter, how¬ 
ever, exists between the halos and, beyond particle mass, 
does not have any properties to call its own. 


Mvir = ^H®(2) RL 


( 1 ) 


where G is the gravitational constant and H( 2 ) is the Hubble 
parameter (Groton et al. 20061. Rvu is the radius beyond 


which the halo’s density falls below 200 x the critical density 
of the Universe. The quoted values of Mvir in this paper 
are calculated by measuring the spherical overdensity within 
the virial radius, rather than the number of bound particles, 
although these values are very similar for all but the smallest 
halos. We refer the reader to the SUBFIND documentation 
for more information. 
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2.2 Considering the Unbound Matter 

Because unbound matter is, by definition, not part of a 
larger object, we consider unbound dark matter particles 
in relation to the halos they are near. Around each halo we 
define a window of radius and count the number of un¬ 
bound particles that fall within that window. If a particle 
falls within more than one window, it is counted as being 
in each window. We do this because we are interested in 
the bulk effect of the surrounding unbound matter and its 
relation to other factors, which are discussed in Section 5. 
We also consider doubly counted particles in Section but 
individual particles are less important for the questions con¬ 
sidered here than the net amount. 

is defined in two different ways throughout this pa¬ 
per. The first method is halo dependent, where we take r^^ 
to be a multiple of the virial radius, Rvir- With this metric, 
the windows around more massive halos will probe larger 
volumes. Since a halo’s gravitational sphere of influence is 
proportional to its radius, a more massive halo influences a 
larger surrounding region, making this a reasonable means of 
defining the window. In the mMR, Rvir is defined in physical 
coordinates and our particles’ locations are in comoving co¬ 
ordinates. To account for this, we translate everything into 
comoving coordinates. 

In the second method, we take r„ to be a fixed value, 
either 1, 2, or 5 h~^ Mpc, so that each halo probes the 
exact same volume regardless of mass. This window is more 
easily comparable to observations than the halo dependent 
window, although it has fewer physical implications for the 
halo. The window radius is in comoving coordinates. 

Figure shows an example of the evolution of bound 
and unbound matter relating to a group sized (10^® h~^ 
Mq) halo from z = 2 until z = 0 (left to right) in a 4 x 
4x4 (comoving) h~^ Mpc window. The green shading and 
contours are the bound matter, defined by particles tagged 
as bound in the halo finder. The purple shading is unbound 
matter, defined as particles not bound to a halo. The solid 
black circle is 1 Rvir, considered the extent of the halo, and 
the dotted circle is 5 Rvir. 

Although the bound matter extends beyond the virial 
radius, in our analysis we assume it fills the volume out to 
1 Rvir and no further. Additional bound matter can be seen 
in the form of smaller nearby halos. Even though they are 
outside the virial radius of the central halo, some of these 
might already be satellites. 

As the halo evolves (left to right) and grows, the un¬ 
bound matter appears to cluster around it, to the point 
where it almost forms a ring circumscribing the bound ma¬ 
terial. We will look at this evolution more closely in Section 

m 

2.3 The Semi-Analytic Model 

To include galaxy properties in our analysis, we post pro¬ 
cess the dark matter halos with a semi-analytic model. We 
ic Galaxy Evolution (SAGE) code from 
, with updates as described in Croton 
et al. (in prep). SAGE models the movements of baryons 
through reservoirs of stars, cold gas, hot gas, and ejected 
gas. Metals are created as part of the life cycle of stars and 
follow the gas, increasing the cooling. The model also in- 


use the Semi-Analyt 

p006 


Croton et al. 


eludes feedback from supernovae and active galactic nuclei, 
which depletes the cold gas reservoir and hinders the cre¬ 
ation of new stars. 

The resolution of the semi-analytic model, when cou¬ 
pled with the mMR, is a stellar mass of ~ 10®'® h~^ Mq. 
To avoid issues relating to this, we use a lower limit of stel¬ 
lar mass of 10® h~^ M©. This creates a sample of 8,411 
galaxies at z = 2 and 9,457 at z = 0. The full sample is 
our density defining population, which we use to measure 
the galaxy count environment. For our analysis, we further 
limit the sample of halos we discuss here to those contain¬ 
ing central galaxies, eliminating satellites, and bringing the 
two populations to 7,089 and 6,870 for z — 2 and z = 0, 
respectively. Satellites are included as part of our density 
defining population when we discuss environment in Sec¬ 
tion 5.1, but are not included in any other capacity, as they 
occupy subhalos, and are therefore unlikely to accrete any 
additional matter. As shown in Croton et al. (20061, SAGE 
agrees well with observations of the stellar mass-metallicity 
relation (e.g Tremonti et al.||2004 1 and the stellar mass-gas 
fraction relation (e.g Cooper et al.|2008l, among others. 


3 UNBOUND MATTER AROUND HALOS 

Halos are the natural result of matter clustering, clumping, 
and virializing. After a halo has formed it can continue to 
grow through mergers and by accreting unbound matter that 
has clustered nearby, pulled in by gravity. When looking at 
the distribution of unbound matter in the universe, it is 
therefore natural to use halos to define the windows where 
we look for unbound matter. We centre a spherical window 
of a given size on each halo. As described above, we measure 
the unbound mass using two different windows, one that 
scales with virial radius, and one that is fixed in size. 


3.1 Halo Dependent Windows 

In Figure we show the virial mass of a halo against the 
amount of unbound mass within volumes defined by radii 
of 2 Rvir, 3 Rvir, and 5 Rvir (left, middle, and right, re¬ 
spectively) at redshifts of 2 = 2 (top) and 2 = 0 (bottom). 
Each point is a halo. The solid line is a fit to the data, with 
a and /3 being the slope and the intersect in the equation 
log^o(Munbound) = « logio(Mvir) + p. The dotted line is the 
mass that would be expected within the window if the un¬ 
bound matter were evenly distributed throughout the box. 

At both redshifts and with all three window sizes, the 
amount of nearby unbound matter closely tracks the virial 
mass of the halo. This is partially due to the window size be¬ 
ing dependent on the central halo, although if the unbound 
matter were randomly distributed, there would be far less 
unbound matter within the window, as evidenced by the 
dotted line. The scatter is very low for all windows. 

A slope greater than one indicates that there is more 
unbound matter around massive halos, relative to their size, 
compared to the regions around less massive halos. A slope 
less than one indicates that there is more unbound matter 
around less massive halos, relative to their size. All of our 
windows result in slopes very close to one, meaning there is 
only a slight trend with halo mass on any scale. 

At both redshifts, the slope of the distribution decreases 
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loSlo(Punbound/p) loglo(Pbound/p) loglo(Punbound/p) loSlo(Pbound/p) 

Figure 1. The evolution of a group sized (10^^ h~^ Mq) halo from z = 2 to z = 0. The green shading and contours represent bound 
matter and the purple shading represents unbound matter. Each panel is a 4 x 4 x 4 h~^ Mpc comoving box. The solid black circle is 
the virial radius of the halo, translated into comoving coordinates. The dotted black line is 5 Rvir- 




logio(Mvir/^ ^Mq) logio(M^i,//i ^Mq) logio(My/i ^Mq) 

Figure 2. Halo mass versus the mass of the unbound matter contained within windows of r^ = 2 Rvir 3 Rvir (middle), and 5 

Rvir (right) centred on each halo. The top row is 2 : = 2 and the bottom row is 2 : = 0. The solid black lines are a fit to the data with the 
parameters of the fit listed in the legend. The dotted black line is the expected mass of the unbound matter within the window if the 
unbound matter were evenly distributed. 
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with increasing radius while the intercept increases. A de¬ 
creasing slope could result from the fact that halos of all 
masses are likely to be within a few virial radii of massive 
halos, as measured by the massive halo, so although there is 
more matter overall near massive halos within that window, 
much of it is bound in other halos. The increasing intercept 
is due to the overall increase in the clustering of matter, 
despite there being fewer unbound particles at 2 ; = 0. The 
closer the fitted line (solid) is to the expected mass (dotted), 
the closer the unbound matter is to the average density. 

From 2 = 2 to 2 = 0, the slope increases at each ra¬ 
dius, although far more drastically (20%) for 2 Rvir than 
for 5 Rvir (1%). This is due to the increase in clustering of 
unbound matter around halos, which correlates to the in¬ 
crease of clustering between galaxies. The change in slope 
with time shows that massive halos pull in disproportionally 
more matter than less massive halos. 

3.2 Fixed Scale Windows 

In Figure]^ we again show the virial mass of the halo against 
the amount of nearby unbound matter, but now we are mea¬ 
suring the unbound matter using a window of fixed radius. 
This method is easier to compare to observations than a 
halo-based window. Left to right, r^ = 1, 2, and 5 h~^ Mpc 
(purple, red, and yellow, respectively) at 2 = 2 (top) and 
2 = 0 (bottom). The solid line is the fit to the data, using 
the parameters a and (3 as described in Section [3.1[ The dot¬ 
ted line is the amount of unbound matter expected within 
the window assuming a uniform distribution. Since the win¬ 
dow is a fixed size, the amount of unbound matter expected 
within it is the same, regardless of halo size. 

At both redshifts and in all three window sizes, there is 
some dependence on halo mass. A slope of zero would signal 
no dependence, but each window has a positive slope, indi¬ 
cating that there is more unbound matter around massive 
halos than around less massive halos when measured using 
the same volume. This is related to the well known relation 


between halo mass and local density (e.g. 

Lemson & Kauff- 

mann|19991|Macci6 et al.|2007||Crain et al.|2009||Haas et al. 

2012 

). The halo mass-environment relation is dependent on 


an environment that is defined by the number of galaxies 
nearby. This figure shows that massive halos not only have 
more halos nearby, they also have more unbound matter. 


There is a halo mass dependent lower limit to the 
amount of unbound matter within a fixed window, which in¬ 
creases with increasing halo mass. It is most striking within 
a 1 h~^ Mpc window, although it is visible out to a radius of 
5 h~^ Mpc. There is a turn-around at Mvu ~ Mq 

in the 1 Mpc window, where the amount of unbound 
matter within the window starts to decrease for cluster sized 
halos. Halos of this mass occupy a large fraction of the I h~^ 
Mpc window, leaving less volume to be occupied by unbound 
matter. 

The strongest correlation between Munbound and Mvir 
comes at small scales, with the 1 h~^ Mpc window having 
twice the slope of the 2 h~^ Mpc window. The slopes are 
fairly stable across redshifts, as are the intercepts. The lat¬ 
ter is an interesting result, since there is roughly half the 
unbound matter at 2 = 0 than at 2 = 2, as shown by the 
decrease in the dotted line for a given window across red¬ 
shifts. It would therefore stand to reason that the intercept 


logio(Mvir) 
in h~^ Mq 

Rvir (2 = 0) 

in h~^ Mpc 

Rvir (2 = 2) 
physical 
coordinates 

Rvir (2 = 2) 
comoving 
coordinates 

10.5 

0.051 

0.025 

0.075 

11.5 

0.111 

0.053 

0.159 

12.5 

0.239 

0.115 

0.345 

13.5 

0.514 

0.247 

0.741 


Table 1. The virial radii of halos at different redshifts. First 
column: halo mass; second column: virial radius at 2 = 0; third 
column: virial radius at 2 = 2 in physical coordinates; fourth 
column: virial radius at 2 = 2 in comoving coordinates. 

should decrease with time. The only scale this is true for is 1 
h~^ Mpc. For the larger scales, the slope decreases slightly 
and the intercept increases, also by a small amount. This 
stability could indicate that the unbound matter is flowing 
into these regions at the same rate at which it is accreting 
on to halos. 

Unlike in the halo dependent windows, there is quite a 
bit of scatter in the relations between Myir and Munbound, 
especially at low redshift. The scatter decreases and an 
upward trending tail emerges for the smallest window at 
masses above 10^^ h~^ Mq at both redshifts, although it is 
much more pronounced at 2 = 0. In the 2 h~^ Mpc win¬ 
dow, the scatter decreases and the tail emerges at a higher 
mass (~ 10^^ ® h~^ Mq). In the largest window, the scat¬ 
ter decreases with increasing Mvir, but an the tail does not 
emerge. A larger simulation with a higher quantity of mas¬ 
sive halos might show the tail for the largest window, but 
the mMR has too few cluster sized halos. 

For comparison purposes between the two window scal¬ 
ings, we have listed the virial radii for several masses at 
both redshifts of interest. They can be found in TableFor 
a 10^^ ® h~^ Mq halo, a 1 h~^ Mpc window corresponds to 
a few virial radii at both redshifts of interest. 


4 FREE-FALL TIMES 

Finding the amount of unbound matter around a halo will 
give some indication of how much that halo will grow by ac¬ 
cretion in the future. Proximity, however, only implies that 
the unbound matter will fall onto the halo on a short enough 
timescale to make a difference in the evolution of the galaxy. 
To measure that timescale, we use the free-fall time of the 
halo plus the nearby unbound matter within our above con¬ 
sidered windows. 

The free-fall time of an object is the amount of time 
it takes for that object to collapse under its own gravity, 
assuming there are no pressure forces countering it. For a 
spherically symmetric object of density p, the free-fall time 
is 

where p = M/(|7r R®). 

For a halo in the Millennium simulation, virial mass and 
radius are related by Equation so the average density of 
a halo is 
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logio(Mvir/^ ^Mq) ^Mg) logio(M^i,//i ^Mg) 



lOgl0(Mvir/^ ^Mg) 10glo(M^i//l ^Mg) 10glo(M^i//l ^Mg) 


Figure 3. Halo mass versus the mass of the unbound matter in windows of 1 (left), 2 (middle), and 5 (right) h~^ Mpc, centred on each 
halo. The top row is 2 : = 2 and the bottom row is 2 : = 0. The solid black line is a fit to the data, with the parameters of the fit being listed 
in the legend. The dotted line is the expected mass of the unbound matter within the window, assuming it is uniformly distributed. 


Mvir 


3 100 
47r G 


n{zf. 


(3) 


Combining Equations and the free-fall time of a halo is 


where Ih is the Hubble time, or the age of the Universe at 
a given redshift, assuming expansion has been constant. 

A reasonable measure of how likely a halo is to grow 
from accreting unbound matter can be found in the free- 
fall time of the halo plus the unbound matter. To calculate 
this, we assume spherical symmetry in the distribution of 
Munbound arouud the halo. While this is clearly not accurate, 
it is our version of assuming a spherical cow. 

We are now considering both the mass within the halo 
and the unbound mass out to n x Rvir, so 


M — Mvir -|- Munbound — Mvir(l -b Munbound/Mvir), (5) 


R = n X Rvir ( 6 ) 

Incorporating these values into p, 

Mvir(l Munbound/Mvir) 1 Munbound/Mvir 

(|™3RL) 

(7) 

and the free-fall time from Equations and is now 


tff = iff 

This relation shows that the free-fall time of the un¬ 
bound matter relies on both the spatial extent of the win¬ 
dow, as well as the amount of unbound matter within that 
volume. At a given radius around a halo of a given mass, a 
larger amount of unbound matter will decrease the free-fall 
time. 

Figure]^ shows the free-fall time of the unbound matter 
enclosed within 2 Rvir, 3 Rvir, and 5 Rvir (left to right) at 
z = 2 and 2 = 0 (top and bottom, respectively) as a fraction 


1 + Munbound/M 

V 


( 8 ) 
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10 11 12 13 14 10 11 12 13 14 10 11 12 13 14 


logio(Mvir/^ ^Mq) ^Mq) logio(My/i ^Mq) 



10 11 12 13 14 10 11 12 13 14 10 11 12 13 14 


logio(Mvir/^ ^Mq) logio(M^i//i ^Mq) logio(M^i//i ^Mq) 

Figure 4. The halo mass versus the free fall time of the halo plus the unbound matter as measured in windows of r^ = 2 Rvir (left), 3 
Rvir (middle), and 5 Rvir (right) at z = 2 (top row) and z = 0 (bottom row). The dotted line is the Hubble time. 


of the Hubble time at that redshift. The dotted line in each 
panel is 1 Hubble time. 

At 2 = 2, all of the unbound matter out to 5 Rvir around 
all of the halos in the sample can accrete onto the galaxies 
within a Hubble time. At z = 0, the unbound matter is 
more clustered, meaning the distribution is steeper, so the 
difference in free-fall times between window sizes is more 
pronounced. All of the galaxies have the potential to accrete 
all of the nearby unbound matter out to 2 and 3 Rvir, but 
only 80% of the galaxies can accrete all of the unbound 
matter out to 5 Rvir- 

As in Section we continue by comparing the halo de¬ 
pendent window to the more observationally motivated fixed 
radius window. To calculate the free-fall time of unbound 
matter within a fixed window, we use Equation and define 
n as 


Tl — rw/Rvir. (9) 

Figure|^shows the free-fall times of the unbound matter 
within fixed windows of radii 1, 2, and 5 h~^ Mpc (left to 
right) at z = 2 and 0 (top and bottom, respectively). The 


dotted line is 1 Hubble time. Note the significant change in 
the range of the y-axis. 

At z = 2, all of the galaxies can accrete all of the un¬ 
bound matter out to 5 h~^ Mpc within a Hubble time. At 
z = 0, however, very few halos have the potential to ac¬ 
crete all of their nearby unbound matter within a Hubble 
time. Only 12% of galaxies can accrete everything within 1 
h~^ Mpc, 1% of halos can accrete everything within 2 h~^ 
Mpc, and no galaxies will accrete everything within 5 h~^ 
Mpc. Since the matter in the 5 h~^ Mpc window includes 
the matter in the 1 and 2 h~^ Mpc windows, some halos can 
still accrete some of the unbound matter within this volume. 
Our calculations reflect the potential of a halo accreting all 
of the unbound matter within the window. 

While only a small percentage of sub-L* halos have the 
potential to accrete all of the unbound matter within 1 h~^ 
Mpc, all of the halos over ~ 10^^'^ h~^ Mq can do so. That 
this turn-over occurs at L* is most likely a numerical coinci¬ 
dence, but it is one we can potentially exploit future analysis. 
A fixed window of 1 h~^ Mpc is a better predictor of future 
growth for galaxies in halos more massive than than 

it is for less massive halos. It is also a better predictor when 
compared to windows of larger radii. 
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Figure 5. The halo mass versus the free fall time of the halo plus the unbound matter as measured in windows of r^ = 1 (left), 2 
(middle), and 5 (right) Mpc. The dotted line is the Hubble time. Note the change in range on the y-axis between 2 = 2 (top row) 
and 2 = 0 (bottom row). 


As mentioned above, in our current model, individual 
unbound dark matter particles can be counted more than 
once if they are near more than one halo. The fraction of 
doubly (or more) counted particles within 5 Rvir of a halo 
is 39% at 2 = 2 and 46% at 2 = 0, with each particle 
found in a window being counted an average of ~ 1.7 times. 
This fraction is even higher for particles within 5 h~^ Mpc 
of a halo, with the fractions being close to 100% at both 
redshifts, and each particle being counted ~37 times at 2 = 2 
and ~70 times at 2 = 0. If each dark matter particle were 
counted only once based on the gravitational pull of each 
nearby halo, the amount of dark matter around a given halo 
would decrease (dramatically for matter calculated within a 
fixed window) and the free-fall times would increase. This 
would affect sub-L* halos more than massive halos, since the 
former have smaller gravitational potential wells and would 
therefore be less likely to pull the particles in. 

The free-fall time only measures the potential of a 
halo to accrete nearby unbound matter and inherently does 
not include any pressure forces or angular momentum that 
would prevent collapse. We are not interested in the col¬ 
lapse of the system, but rather the accretion of the unbound 


dark matter particles onto the halo, so the free-fall time is a 
proxy for a more complex measurement which considers the 
time it would take for a given particle to reach the edge of 
a halo. We expect this to correlate with the infall of diffuse 
gas that could ultimately go on to fuel future star formation 
and galaxy growth, and hence we can speculate about the 
fate of the galaxy itself. 


5 OBSERVABLES 

The correlation between halo mass and the amount of nearby 
unbound matter measured in Section|^shows that halo mass 
is a reasonable predictor of nearby unbound matter, depend¬ 
ing on window type and scale. Unfortunately, halo mass is 
not an easily measurable quantity observationally. This, in 
turn, makes Rvir an unreliable basis for a metric. To make 
our results more comparable to observations, we now switch 
to using only a fixed window for measuring the unbound 
matter and extend the previous analysis to compare the 
nearby unbound matter with mock galaxy properties. From 
this we can make predictions of how diffuse matter is dis- 
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tributed which can ultimately be tested by observers using 


instruments, such as the Cosmic Web Imager (Martin et al. 


2014), or new telescopes, such as the SKA, to measure the 


intergalactic gas which traces the unbound matter. 


5.1 Galaxy Environment 


An obvious first “observable” to look at is the environment 
of the galaxy as measured by the nearby galaxies. We use a 
spherical fixed aperture metric to define environment, which 


has been discussed extensively in Muldrew et al. (2012). To 
avoid confusion, we will continue to use the word “window” 
to describe the volume containing unbound matter and use 
“aperture” to describe the volume over which we are count¬ 
ing galaxies. It is useful to note that windows are centred 
on the halos and the apertures are centred on the galaxies. 
Within a semi-analytic model, these are the same location, 
but this is not always true in hydrodynamic simulations or 
in observations. 

A galaxy or halo’s local density is defined as 


p/p = 


N, 


4/37rr5^ p 


( 10 ) 


where Ng is the number of galaxies within the aperture (in¬ 
cluding itself, other central galaxies, and any satellite galax¬ 
ies that fall within the aperture), r^ is the radius of the aper¬ 
ture, and p is the mean number of density defining galaxies 
per unit volume. 

For reference, a single galaxy in our simulation within a 
1 h~^ Mpc aperture has a p/p = 6.16, in a 2 h~^ Mpc aper¬ 
ture p/p — 0.77, and in a 5 h~^ Mpc aperture p/p — 0.05. 
By definition, all galaxies are over-dense wh en measured on 
a 1 h~^ Mpc scale. As has been discussed in 


120121; Haas et al. (2012 I; Shattow et al. (20131, etc., the en¬ 


Muldrew et al. 


vironment of a galaxy, as measured by galaxy counts within 
an aperture, correlates with the mass of the halo containing 
the galaxy. 

Figure shows the fixed aperture galaxy count of the 
galaxies against the amount of unbound matter within aper¬ 
tures of r^ = 1, 2, and 5 h~^ Mpc, shown by the purple, 
red, and yellow dots, respectively. Galaxy-count environ¬ 
ment is by nature a quantised measure and we do not make 
allowances for part of a galaxy being within an aperture, 
resulting in discrete values for the density contrast at low 
density. These two quantities are closely correlated on all 
three scales. The largest radius aperture/window has the 
tightest correlation and the smallest radius has the weakest. 

To put the environment versus Munbound into the same 
context as our previous figures, we show in Figure the 
velocity dispersion of the halo against the mass of unbound 
matter within windows of r^ = 1, 2, and 5 h~^ Mpc (left to 
right) at 2 = 0. We shift from using the halo’s mass to its 
velocity dispersion, a, which is the 1-dimensional velocity 
dispersion of the halo particles. It is related to the velocity 
dispersion of galaxies in a group or cluster and is therefore 
a more observable property. The black line is the same fit 
to the data as in Figure with the substitution of a for 
Mvir, and the dotted line is the expected amount of unbound 
mass in that window. The colour is based on the log of the 
density contrast, measured in apertures the same size as the 
windows. 
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Figure 6. The local density contrast, as measured by the number 
of galaxies within fixed apertures of 1 (purple), 2 (red), and 5 (yel¬ 
low) h~^ Mpc, versus the mass of the unbound matter contained 
within windows of the same radius at 2 = 0. 


The amount of unbound matter depends strongly on the 
environment of the halo and only weakly, if at all, on the halo 
velocity dispersion. The correlation between the halo and the 
nearby unbound matter measured in a fixed window appears 
to result largely from the correlation between halo mass and 
environment. At 2 = 2 the same pattern emerges but with 
less scatter, similar to the difference in redshifts of Figure 
The halos also inhabit a smaller range of environments at 
higher redshift, which is predicted from structure formation. 

We see the same pattern in Figure which is the same 
as the bottom panel of Figure]^ with a instead of Mvir on 
the x-axis and is now coloured by the galaxy count environ¬ 
ment as measured in apertures with radii the same size as 
the unbound matter windows, namely 1, 2, and 5 h~^ Mpc 
(left to right, respectively). Here, the potential future evo¬ 
lution of the halo is highly dependent on the environment, 
with lower density halos (blue points) requiring far more 
time to accrete the nearby unbound matter as compared to 
halos in high density environments (red points). This comes 
directly from the calculation of tn in Equation]^ where we 
showed the free fall time is dependent on both the amount 
of nearby unbound matter and the distance out to which the 
unbound matter extends. 


5.2 Stellar Mass and Other Galaxy Properties 

Thus far we have limited our measurements to the virial 
mass and velocity dispersion of a galaxy’s halo, i.e. galaxy 
properties that emerge from the dark matter simulations 
alone. We now include observable properties from the bary- 
onic physics in the semi-analytic model to test their inter¬ 
dependence on the amount of nearby unbound matter. We 
also use unbound matter as a proxy for the diffuse gas of 
the IGM, which can potentially be observed. To this end, 
we switch the y-axis from Munbound to ft x Munbound to rep¬ 
resent the diffuse gas, in a similar manner to [Shattow et al.| 
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Figure 7. Velocity dispersion versus the mass of the unbound matter within fixed windows of 1 (left), 2 (middle), and 3 (right) h~^ 
Mpc at z = 0. These are the same distributions as the bottom row in Figure]^ with velocity dispersion rather than Mvir- The colour is 
the local density contrast, as measured by the number of galaxies within a fixed aperture of r^ = 1 (left), 2 (middle), and 3 (right) h~^ 
Mpc. The solid black line is the fit to the data and the dotted black line is the expected mass of the unbound matter if it were evenly 
distributed. 



logio((T/km s ^) logig((j/km s ^) logjg(cr/km s ^) 

Figure 8. Velocity dispersion versus the free fall time of the halo plus the unbound matter within fixed windows of 1 (left), 2 (middle), 
and 3 (right) h~^ Mpc at 2 = 0. These are the same distributions as the bottom row in Figure]^ with velocity dispersion rather than 
Mvir- The colour is the local density contrast, as measured by the number of galaxies within a fixed aperture of r^^ = 1 (left), 2 (middle), 
and 3 (right) h~^ Mpc. The dotted line is 1 Hubble time, or the age of the universe at 2 = 0. 


The baryon fraction, is assumed to be the cos¬ 
mic value of 0.18 ( Planck Collaboration et al.||2015 l. Semi¬ 
analytics show a baryon fraction for galaxies to be slightly 
higher than this (Croton et ah, in prep), leaving the IGM a 
lower fraction overall. It is likely, however, that the deficit 
is found in regions of very low density, rather than in the 
volumes directly surrounding the halos. This results in the 
cosmic baryon fraction being a good estimate in the regions 
of unbound matter around halos. 


Figure shows the stellar mass of a galaxy against the 
amount of diffuse gas (i.e. the unbound matter multiplied by 
the baryon fraction), within a.2 h~^ Mpc window. The black 
lines are the fit to the data, with the slope and intercept 
listed in the top panel. All three panels are identical except 
for the colour coding, which represents a different observable 
in each panel: (top) the galaxy-count environment measured 


in a 2 h~^ Mpc aperture, (middle) the gas fraction of the 
galaxy, and (bottom) the metallicity of the galaxy disk, all 
at 2 = 0. 

Considering first all three panels together, the diffuse 
gas mass has a lower limit which increases slightly with stel¬ 
lar mass, as is the case in Figure although it is not as 
striking as the halo mass dependent limit there. The slope 
of the distribution is 0.10, ~ 30% less than the slope for 
the comparison against Mvir, using the same window. The 
scatter is also significantly larger with M*. Stellar mass, un¬ 
like environment, is therefore a poor predictor of unbound 
matter within a fixed window. 

The top panel, colour-coded to show galaxy-count en¬ 
vironment, reflects the findings of Figures and [T] above, 
namely that galaxy-count environment closely tracks the 
amount of unbound matter, regardless of galaxy mass. In 
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other words, for a fixed stellar mass, the environment and 
the amount of nearby diffuse gas are strongly correlated. 

The middle is colour-coded to show the gas fraction 
within the galaxy, where fgas = Mgas/(M* + Mgas). We 
have tagged the galaxies with gas fractions of less than 1% 
with pink stars. There is no discernible trend in these galax¬ 
ies. Similar figures tagging the high gas fraction galaxies 
have shown the same result. 

The bottom panel is colour-coded to show the metallic- 
ity of the cold gas in the disk of each galaxy. Metallicity is de- 
hned as Z = Mz/McoidCaa and is in units of solar metallicity, 
Zq= 0.02 i Grevesse fc Sauval[]1998 |. Low metallicity galax¬ 
ies are tagged with stars and show no trend with Munbound- 
The excess of low metallicity galaxies at low masses is due 
to the fact that new galaxies in the model form out of pris¬ 
tine gas and have not had much time to create metals. Like 
in the middle panel, there is no dependence on the enclosed 
diffuse gas mass. 
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Figure 9. Stellar mass versus the mass of the unbound baryonic 
matter found in fixed windows of radius r.^^ = 2 h~^ Mpc at 
z = 0. Top panel: colours are the local density contrast measured 
by the number of galaxies within a 2 h~^ Mpc fixed aperture. 
Middle panel: colours are the galaxy gas fraction. Galaxies with 
fgas < 0.01 are tagged with stars. Bottom panel: colours are the 
metallicity of the galaxy disk. Galaxies with [Z/Zq] < —1 are 
marked with stars. 


6 DISCUSSION 
6.1 Unbound Matter 

At the beginning of this paper, we asked three questions 
regarding the unbound matter in the universe. The first, 
regarding its distribution and evolution over time, is an¬ 
swered in Figures[^-[^ When measured in a halo dependent 
window, the relation between the virial mass of a halo and 
the amount of unbound matter nearby becomes steeper at 
low redshift, indicating that unbound matter has clustered 
around the more massive halos at a faster rate than around 
less massive halos, confirming previous studies. 

The Mvir-Munbound relation is also slightly steeper for 
smaller windows, i.e. r.^^ = 2 Rvir, than it is for larger win¬ 
dows, i.e. r.^,y = 5 Rvir, with the difference being more pro¬ 
nounced at 2 = 0 than at higher redshift. This narrowing of 
the distribution mimics the growth of bound structure, with 
the halo filling the nearby volume with material pulled in 
from the outer volume. There is a slight increase in scatter 
in the relations, both with increasing radius and with time, 
with the largest amount of scatter being measured at 5 Rvir 
from the halo at 2 = 0. 

The dark matter directly outside of the halo is con¬ 
trolled by two opposing forces - the gravity pulling it to¬ 
wards the halo, and the galaxy feedback processes pushing 
it away. Since our model is built on an N-body simulation 
and does not include feedback processes until we add the 
semi-analytic model, we are limited to analyzing the effects 
of gravity alone. Most of the literature on the interplay be¬ 
tween baryons and halos focuses on the inner part of the 
halos, where baryons have the strongest effect and their in¬ 
fluence is more easily observed. |Duffy et al.| ( [^10[ ), however, 
show the outer edges of galaxy sized halos are slightly af¬ 
fected by feedback, with weak supernova feedback decreas¬ 
ing the amount of dark matter at the edge and strong feed¬ 
back from AGN or supernovae increasing it. It is unlikely 
that these feedback processes would push the unbound mat¬ 
ter more than an additional virial radius beyond the halo, 
even in a small halo. The amount of unbound matter, includ¬ 
ing both ejected and not-yet-bound material, would there¬ 
fore increase in the annulus between 1 and 2 Rvir, but not 
affect the rest of the distribution by much. 
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Unlike in the halo dependent window, the slope of the 
relation between unbound and virial mass does not change 
much with redshift for a given fixed-size window. The scat¬ 
ter, however, does increase considerably. Both the scatter 
and the slope of this relation decrease with increasing ra¬ 
dius, as the volume probed by the fixed window approaches 
the mean density. The scatter is especially prevalent aronnd 
smaller mass halos, reflecting the variety of environments 
that low mass halos inhabit. For example, the group sized 
halo (~ 10^® h~^ Mq) shown in Figurej^has several dwarf 
galaxy halos and subhalos within a few virial radii. This 
large halo does not live within the spheres of influence of 
the low mass halos, but it does live within a few Mpc. 
These dwarf galaxy sized halos, therefore, will have a similar 
amount of unbound matter around them as the group sized 
halo. 

The virial radius is defined as the point at which a den¬ 
sity fluctuation is 200 times the critical density. This is a 
somewhat arbitrary definition (see Shull 2014 and refer¬ 
ences therein for a full discussion), since the slope of the mass 
profile continues outside of the determined virial limit, as 
shown in Figure T| If th e virial radius were determined to be 
smaller, as Shull (20141 finds for the Milky Way, the amount 
of matter that would be considered unbound in these sim¬ 
ulations would increase. If Rvir were defined as the point 
where a halo reaches 180 x pcrit, another common defini¬ 
tion, its value would be larger and the amount of unbound 
matter would be lower. Likewise, if the boundary of 200 
times the mean mass density (flm) is used, the amount of 
unbound matter decreases even further, as discussed by 


Daalen & Schaye (20151. Also, as discussed in Section 


2.1 


the amount of unbound matter changes depending on which 
halo finder is used. These profiles and relations are therefore 
dependent on several different factors, which are more likely 
to affect the results of the halo based window than the fixed 
radius window. 


6.2 Free-fall Times 

The second question we ask, regarding the likelihood of a 
halo to accrete all of the unbound matter within a window, 
we answer using the free-fall times, as shown in Figures|^and 
0 The free-fall time is calculated for the matter in the halo 
plus the unbound matter within a certain window. Lower 
density regions have larger free-fall times, and the farther 
from the halo a window reaches, the lower the total density 
will be. This accounts for the increase of free-fall times for 
larger windows. The increase in scatter with window size 
and decreasing redshift is magnified in the free-fall times. 

Other effects to consider include the fact that halos have 
grown between z = 2 and 2 = 0, so the halo dependent 
window around each halo covers more volume. To counter 
this, we are using comoving coordinates, which decrease the 
extent of the radius with redshift by a factor of (1 -|- z)~^ 
relative to the physical extent, as shown in Tablej^ The solid 
black circle representing 1 Rvir in Figurej^shows the growth 
of the halo is the dominant effect of these two between 2 = 2 
and 2 = 0.5, although the radius decreases again between 
2 = 0.5 and 2 = 0. 

We have marked the Hubble time on Figures and 
showing what fraction of halos will be able to accrete all 
of the unbound matter within a given window in a Hubble 


time, as measured at that redshift. It is also interesting to 
consider how many of the halos at 2 = 2 could have accreted 
all of their enclosed unbound matter by 2 = 0. tH(z = 0) ~ 
3 X tH(z = 2), so the unbound matter would have to be 
within 2 additional Hubble times (as measured at 2 = 2) 
for the halo to be able to accrete all of the enclosed matter. 
From these figures, we can see that the free fall times of the 
enclosed masses out to both 5 Rvir and 5 h~^ Mpc are less 
than 1 fn at 2 = 2, well within the cutoff for accretion. 

As mentioned above, the free-fall time measures the 
pressure-less collapse of a spherical distribution. Since we 
do not have a pressure-less system and it is obvious from 
Figure that we do not have a spherical object, it is in¬ 
teresting to note how many of these particles are actually 
accreted by halos between 2 = 2 and 2 = 0. 

By comparing the particles that are within 5 Rvir of 
a halo at 2 = 2 to the particles that are part of a halo 
at 2 = 0, we find that 78% of nearby particles are accreted. 
The same comparison to particles within 5 h~^ Mpc at 2 = 2 
shows that only 36% of these proximal particles are accreted 
in the intervening time. The lower accretion rate using the 
fixed window size most likely results from the large number 
of multiply counted unbound particles in that metric, which 
increases the free-fall times when the particles are counted 
only once (see Section . Future work will take this into 
account for a more rigorous calculation of the free-fall times 
for the purpose of predicting future evolution of a halo. 

If we assume baryons trace the dark matter particles as 
halos accrete them, then the very long free-fall times (which 
are likely lower limits on the accretion time-scales) that have 
emerged by 2 = 0 could indicate a natural quenching of 
the growth of the interior galaxy that has nothing to do 
with baryonic processes. Halos are merely accreting mate¬ 
rial more slowly, which will increasingly slow down the star 
formation cycle on which the galaxy depends. In this case, 
semi-analytic models might assume too high an efficiency 
for feedback processes when, in reality, halos might quench 
themselves. This idea is partially backed up by Figure 
which shows that the local galaxy count density might play 
a large role in the future evolution of the galaxy. 


6.3 Observables 

Our third question addresses if there are any correlations 
between observable properties of a galaxy and the amount 
of unbound matter, with the intention of using these corre¬ 
lations to make predictions for the future observation of the 
unbound structure of the Universe. 

We first look at the unbound matter against the galaxy 
count environment, both on its own and as a function of 
the velocity dispersion of the halo. These two properties are 
closely related, with the strongest correlation at a scale of 
I'w ~ ~ ^ Mpc. Thus, a 5 h~^ Mpc aperture is a 

good means of determining how much unbound matter is 
inside a window of the same size, regardless of galaxy mass. 
Taking Figure|^into account, however, shows that this scale 
probes very close to the mean density of unbound matter. 
Its usefulness as a probe of future growth is therefore less 
certain. 

From Figure we see that zero halos at 2 = 0 have 
the potential to accrete all of the unbound matter out to 
5 h~^ Mpc by the time the universe doubles in age, i.e. 
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within another Hubble time as measured at z = 0. Even 
on a 1 h~^ Mpc scale, only 12% of the galaxies can accrete 
all of the unbound matter on this timescale, making the 
smaller radius a more physical choice of window/aperture 
size to follow to determine future growth, although having 
a weaker correlation and more scatter means it is less clear 
how much unbound matter surrounds the halo based on its 
galaxy-count environment. 

Interestingly, halos of all masses with more nearby 
galaxies also have more unbound matter, and therefore 
shorter free fall times. Rather than depleting the unbound 
matter by accreting it all onto bound halos, for a halo of a 
given mass having more nearby galaxies is a strong indicator 
of relatively more unbound mass within the window. This is 
especially apparent in Figures |^-[^ 

This effect is even more pronounced when galaxy stellar 
mass is compared to the diffuse IGM in Figure There is 
little correlation between stellar mass and diffuse gas in a 
fixed radius window, since galaxies of a given mass can oc¬ 
cupy a range of halo masses, even as the central galaxy. Any 
pattern in the halo mass and the unbound dark matter from 
Figureis therefore washed out. The galaxy count environ¬ 
ment, however, shows a clear trend with unbound matter in 
the top panel of[^ Of all of the possible observable properties 
we can measure using the model, the galaxy count environ¬ 
ment is the best indicator of how much unbound matter is 
within a window of the same size. 

Different calculations of environment often correlate 
one another (e.g. [Muldrew et al.|[2012[ [Shattow et al.| 
I because their goal is to probe the underlying structure 
of the mass density field. Most environment metrics are lim¬ 
ited to observables, such as the number of galaxies nearby, 
the distance to the Nth nearest neighbour, or the amount 
of volume closer to a galaxy than to any other galaxy (i.e. 
Voronoi tessellation). These metrics all use galaxies, or the 
peaks of the density field as their density defining popula¬ 
tions. This work shows that the relation between environ¬ 
ment metrics continues to lower density regions of space - 
where the diffuse gas lives. 

Gas fraction and metallicity are also associated with 
stellar mass, and to a lesser extent, environment (Thomas 
et al. I |2005| |De Lucia et al.||2006| [Cooper et al.||2008| |Davq 

et al.|2011[ ). These two galaxy properties, however, show no 

correlation with nearby unbound matter beyond their initial 
dependence on stellar mass. Since they do correlate with 
environment, one might expect some variation with unbound 
matter. The lack of a trend is possibly due to the type of 
galaxy model used. Semi-analytic models do very well at 
reproducing observable properties and relations of galaxies, 
but they are based entirely on the halos. They therefore 
lack environment information beyond the virial radius of the 
halo in which the galaxy resides. A more complex model or 
a hydrodynamic simulation might be required to tease out 
these effects. 

Overall, galaxy count environment is the closest mea¬ 
surement we have to a smoking gun of nearby unbound mat¬ 
ter, especially when they are both measured on the same 
scale. 


7 SUMMARY 

This paper addresses the often-ignored unbound matter in 
simulations in the context of how a galaxy’s (or halo’s) en¬ 
vironment is affected by its evolution. We investigate the 
relationship between the amount of unbound matter and 
halo mass using a variable and a fixed window and examine 
how these relationships evolve over time. We also consider at 
the potential of a halo to accrete the nearby unbound mat¬ 
ter within a Hubble time. Finally, we compare the unbound 
matter to observable properties of galaxies in the hopes of 
finding a clear signal of nearby unbound matter for observers 
to look for. 

The main conclusions of this paper are: 

• Unbound matter clusters around bound matter, with 
more unbound matter being found near massive halos when 
measured on a fixed scale, as seen in Figure For halo 
dependent windows, this is only true from r^ = 2 Rvir at 
z = 0. The remaining windows at both redshifts have slopes 
equal to or less than 1, meaning there is no preference for 
more massive halos. This is shown in Figure 

• When measured in windows scaled to the virial radius 
of the halo, the free-fall time of the total mass inside the 
window reveals that most halos have the potential to accrete 
all of their surrounding unbound matter within a Hubble 
time out to a radius of 5 Rvir- We show this in Figure]^ 

• At z = 2, all galaxies can accrete all of the diffuse mat¬ 
ter out to a fixed radius of 5 h~^ Mpc within a Hubble time. 
At z = 0, only 12% can accrete all of the unbound matter 
out to 1 h~^ Mpc, as seen in Figure 

• Galaxies in denser environments are more likely to have 
more unbound mass nearby, regardless of stellar or halo 
mass, which we show in Figures [^H] and|^ 

• Other galaxy properties, such as gas fraction and disk 
metallicity, show little correlation with the diffuse gas, mak¬ 
ing galaxy-count environment the most promising observa¬ 
tional probe of the large scale distribution of the intergalac- 
tic medium in our analysis, as is clear in Figure]^ 

Future work will look at the evolution histories of spe¬ 
cific halos, measuring how important these accretion events 
are in comparison to other means of growing, namely major 
and minor mergers. This will further address the question 
of how a galaxy’s evolution is affected by its environment. 
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